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We present a novel type of differential equations for on-shell loop integrals. The 
^ ■ equations are second-order and importantly, they reduce the loop level by one, so 

that they can be solved iteratively in the loop order. We present several infinite 
series of integrals satisfying such iterative differential equations. The differential 
operators we use are best written using momentum twistor space. The use of the 
latter was advocated in recent papers discussing loop integrals in Af = 4 super 
Yang-Mills. One of our motivations is to provide a tool for deriving analytical 
results for scattering amplitudes in this theory. We show that the integrals needed 
for planar MHV amplitudes up to two loops can be thought of as deriving from a 
single master topology. The master integral satisfies our differential equations, and 
so do most of the reduced integrals. A consequence of the differential equations is 
that the integrals we discuss are not arbitrarily complicated transcendental func- 
tions. For two specific two-loop integrals we give the full analytic solution. The 
simplicity of the integrals appearing in the scattering amplitudes in planar M = 4 
super Yang-Mills is strongly suggestive of a relation to the conjectured underlying 
integrability of the theory. We expect these differential equations to be relevant 
for all planar MHV and non-MHV amplitudes. We also discuss possible extensions 
of our method to more general classes of integrals. 
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1 Introduction 



In this paper we present new differential equations for on-shell loop integrals. Our main 
motivation is to develop efficient methods to determine the loop-level S-matrix in M = 4 
super Yang-Mills (SYM), which is built from precisely the types of integrals we discuss 
here. The planar M = 4 theory is believed to be governed by some underlying integrable 
structure and hence it certainly deserves to have a beautiful and simple S-matrix. Our 
differential equations can be viewed as a step towards explaining this simplicity, as they 
restrict the integrals contributing to the amplitudes. Our results are at the level of specific 
integrals and so can also be used for subsets of integrals appearing in less supersymmetric 
theories. We also expect our method to be applicable to larger classes of integrals. 

The simplicity of the loop integrals implied by the differential equations we find is sug- 
gestive of a connection to the expected underlying integrability of planar M = 4 SYM, and 
more concretely, of the underlying symmetries of its tree-level amplitudes. It was found 
that the tree-level S-matrix in this theory has a hidden symmetry, dual superconformal 
symmetry [TJ [2]. Together with the ordinary superconformal symmetry of the theory, it 
leads to a Yangian symmetry [3] . At loop level, a priori the symmetries are broken. 

The breaking of the bosonic part of dual superconformal symmetry is well understood 
and is controlled by anomalous Ward identities jH [1] that were initially derived for cer- 
tain Wilson loops, and are conjectured to hold for the scattering amplitudes. The Ward 
identities completely fix the functional form of the four-point and five-point amplitudes, 
confirming the ABDK/BDS conjecture for these cases [5]. Starting from six points, a mod- 
ification is required, called the remainder function [6 J . An alternative way of understanding 
the dual conformal symmetry at loop level is possible by using the massive regulator of [7] , 
which is inspired by the AdS/CFT correspondence [8]. This setup is expected to make dual 
conformal symmetry exact, i.e. unbroken, at loop level. Given this, it is natural to assume 
that the loop-level integral basis of M = 4 SYM should consist of integrals having this 
exact symmetry (see [9] and references therein). Indeed, the absence of triangle sub-graphs 
at one loop was confirmed in [10] and further support for this conjecture comes from [TTj . 

A similar understanding of the full Yangian symmetry at loop level is still lacking and 
subject to ongoing research [12]. Recently it was argued that the (unregulated) loop inte- 
grand has the full Yangian symmetry, up to total derivatives [13J . Indeed the integrand can 
be recursively constructed via a BCFW type recursion relation [13] (see also [E]). Ignoring 
regularisation, this recursive construction can be written as a sequence of Yangian invariant 
operations on the basic Yangian invariant functions [T3] (see also [T5] for a discussion of 
Yangian invariants). 

Just as the dual conformal symmetry of the integrand [16] was a hint that there is an 
anomalous Ward identity [I] , the existence of the Yangian invariant integrand indicates that 
there should be a way to directly understand the breaking of the full symmetry. Although 
we do not yet make a direct connection to the Yangian generators, we find it likely that 
the underlying Yangian symmetry is related to the differential equations we find. 
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Further very concrete motivation for our study also comes from recent explicit results 
for the hexagonal light-like Wilson loop, which is dual to the six-gluon MHV amplitude [B] 
(for reviews see [17] and references therein, and [18] for recent developments). In [19], a 
remarkably simple form of the six-point remainder function [6] was given, based on previous 
work [20]. Simple results are also conjectured to hold in special kinematics [21]. Recently, 
the analytic six-point Wilson loop result of [20] could be reproduced in a kinematical limit 
by an analytic calculation of the corresponding scattering amplitude, extending previous 
calculations for mass-regulated amplitudes [HI 122] • Related recent work on loop amplitudes 
in M = 4 SYM can be found in [23]. See also [24] for other related recent developments. 

The use of differential equations to evaluate loop integrals [25] is widespread in the 
literature, for a review see chapter 7 of [22] ■ In this approach one considers a set of master 
integrals, and an associated family of integrals is be obtained from by shrinking lines. One 
differentiates the master integrals with respect to kinematical invariants or masses. The 
result is in general a linear combination of several integrals. Integral reduction identities are 
then applied to re-express the latter in terms of the set of master integrals. In complicated 
cases this step can be non-trivial, as in general it requires the knowledge of all reduction 
identities. In general one obtains a set of first-order differential equations for the master 
integrals. A disadvantage is that in order to solve for a given integral, one has to deal with 
all integrals of a given family. It can also happen that when computing a finite integral, 
intermediate steps contain divergences that only cancel at the end, and a regulator has to 
be used. 

The differential equations we obtain here are quite different in nature. There are two 
important differences to the method described above: Firstly, our differential equations 
can be applied directly to a given integral, without having to know all integrals with fewer 
propagators, or any integral reduction identities. Secondly, they are second-order equations, 
and, importantly, they reduce the loop level by one. In other words, the homogeneous term 
of the equations corresponds to lower-loop integrals, and the equations therefore have an 
iterative structure. 

While the differential equations method known in the literature can be applied in prin- 
ciple to all loop integrals, the differential equations we find are specific to a certain class 
of integrals. The latter constitute a subset of integrals needed for computing scattering 
amplitudes in a generic theory. There is reason to believe that in the special case of M = 4 
SYM, all planar loop integrals are constrained by the type of differential equations we find 
here. Very concretely we observe that the one- and two-loop MHV amplitudes in Af = 4 
SYM in the form given in |13j . can all be expressed in terms of a single 'master' integral, 
namely 




(1.1) 



at one- and two-loop order, respectively. The reason is that reduced topologies, such as 
box, double box or penta-box topologies, can be obtained by taking soft limits of the master 
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integrals. We will give the precise definition of the integrals later. In the six-point case, 
they were used recently to obtain an analytic result for the remainder function at two loops 
(in a kinematical limit) by two of the present authors in [27]. This is the first time that 
such an analytical result was obtained directly for the scattering amplitude, and the results 
also hinted at a certain simplicity of the integrals that had to be computed. We will make 
this more precise in this paper. As we will show, the integrals depicted in (11.11) satisfy 
second-order differential equations that reduce their loop degree by one, namely 




The one-loop integrals appearing in (II. ip and on the r.h.s. of (II. 2p satisfies a similar 
equation with simple rational functions as inhomogeneous terms. These equations imply 
powerful constraints on the functional dependence of those integrals. In the following we 
will give several examples where we solve such equations and obtain analytic answers for 
two-loop integrals with multiple external legs. 

Iterative differential equations of the type (II. 2p were previously found in [16] by two of 
the present authors and Smirnov and Sokatchev for certain off-shell ladder (i.e. multiple 
box) integrals, and are closely related to the equations we find here. We will review the 
differential equations of [16] in section [3j The differential equations presented here can be 
thought of as a generalisation of the ones used in [16] to on-shell integrals. 

We find it very convenient to work with momentum twistor variables [28J. These are 
well-suited to describe planar loop integrals. They solve the momentum conservation and 
on-shell constraints and therefore are unconstrained variables. This makes them the natural 
variables to use if one is interested in differential equations. One-loop box integrals were 
discussed in momentum twistor space in |29j. In the cases where infrared divergences are 
present we employ the AdS-inspired mass regulator of [7J. The latter allows us to stay in 
four dimensions and continue to use momentum twistors in those cases as well. For recent 
references using momentum twistors and this regulator see [221 EZ] • 

The outline of the paper is as follows. In section [21 we recall the definition of momentum 
twistor variables and the expression for the one- and two-loop MHV amplitudes in M = 4 
SYM. It is shown that at each loop level there is only one 'master' topology. In section [3] we 
show that certain classes of integrals satisfy second-order differential equations that reduce 
the loop order by one. We first review the differential equations for off-shell integrals found 
in [16] , and then generalise them to the on-shell case. We present two mechanisms for finding 
such differential equations. We present several infinite classes of integrals satisfying iterative 
differential equations. In section H] we give an example for how to solve the differential 
equations, using certain assumptions about boundary conditions. We give explicit analytical 
results for several multi-leg integrals at two loops. In section we conclude and comment 
on several possibilities of extending our method. There are two Appendices. In Appendix 
A we discuss twistor differential operators that annihilate the integrand of the integrals 
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appearing in the MHV amplitudes, up to anomalies. Appendix B contains the analytic 
formulas for the one- and two-loop penta-box integrals with "magic" numerator. 

2 Motivation: MHV amplitudes 

In this section we recall how momentum twistor variables can be used to describe loop 
integrals. We also recall the recently-proposed integral representation for planar two-loop 
MHV amplitudes in M = 4 super Yang-Mills |13j . We show that at each loop order, all 
integrals contained in the MHV amplitudes can be thought of as deriving from a single 
master topology. Reduced topologies are obtained by taking soft limits. This is related to 
the consistency of the loop integrand with soft limits [27J. 

Let us briefly recall how the momentum twistor variables are related to the standard 
momentum space variables. Given the n incoming light-like momenta of a planar colour- 
orderd ordered amplitude, 

pf = , (2-3) 
we define the dual coordinates in the usual manner [3U1 US] , 

xf-x%!=pf- (2-4) 

The dual Xi define a light-like polygon in the dual space. On the other hand, a point in dual 
coordinate space corresponds to a (complex, projective) line in momentum twistor space. 
Two dual points are light-like separated if the corresponding lines in momentum twistor 
space intersect at some point in momentum twistor space. Thus the light-like polygon in 
dual space corresponds to a polygon in momentum twistor space with each line intersecting 
its two neighbouring lines as each dual point is light-like separated from its two neighbours. 
The n momentum twistors associated to this configuration of n light-like lines are defined 
via the incidence relations, 

%i — (Ki^i)' — x i a ^ia = Xi+lKa . (2.5) 

The momentum twistor transforms linearly under the action of dual conformal symmetry, as 
indicated by the fundamental s/(4) index A. Moreover the n momentum twistors describing 
the polygon are free variables, in contrast to the dual points Xi which obey the constraints 
of light-like separation from their neighbours. The dual point x^ is associated to the line 
described by the pair (Z^iZi) or (i — 1 i) for short. The incidence relations ( 12. 5 p allow one 
to express functions of the Xj in terms of momentum twistors. For example we have 

2 (j-lij-lj) 

where the four-brackets and two-brackets are defined as follows, 

(ijkl) = e ABCD ZfzfZ^ZP, (ij) = X^Xja . (2.7) 

The four-brackets are obviously dual conformal invariants while the two-brackets are in- 
variant under just the Lorentz and (dual) translation transformations. 
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(2.6) 





Pi-l' 



Zi Zi_i 



Figure 1: 'Two-mass-hard' integral in momentum space, dual space, and momentum twistor 
space variables, see equations (|2.8p and (|2.9p . 



We can consider integrals over the space of lines in momentum twistor space (AB). 
As we have discussed, this is equivalent to an integration over points in dual space. For 
example, the two-mass hard integral of Fig. [1] is given by 



d 4 k (p 3 + . . . + p i _ 1 ) 2 (p 2 +P3) 2 f d^xo x\x\^ 



in 2 k 2 (k + p 2 ) 2 (k + p 2 + p 3 ) 2 (k + p 2 + ■ ■ ■ + Pi-1) 2 J m 2 ^02 x 03 x 04 x 0i 



„2 ' 



(2.81 



where the first expression is written in momentum space and the second one in dual coor- 
dinates. In twistor notation, this becomes 

d^Z AB (23 i— 1 i) (1234) 



m 2 (ABU) ( AB23) {ABM) {AB i-1 i) ' 

The integral ( 12. 8ft as written is infrared divergent, so its proper definition should involve a 
regulator. We will use the the AdS regularisation introduced in [TJ. The latter allows us 
to stay in four dimensions and regulates the infrared divergences by masses that in turn 
are generated by a Higgs mechanism. For actual calculations we will use the regularisation 
where all masses are equal. In particular this means the outermost propagators in the 
planar loop integrals that we are studying are modified as follows, 

\— » 2 I 2 - ( 2 - 10 ) 

X ij X ij + m 

In twistor space this has the effect that each of the propagator factors becomes 

1 1 - 1 

(AB i-1 i) > (ABi-li) = (ABi-li) +m 2 (AB)(i-li) ' 

When doing this, there is a choice of adding 0(m 2 ) terms to the integrand. This is certainly 
relevant if one wishes to obtain the exact m-dependence of amplitudes on the Coulomb 
branch of M = 4 SYM. In the cases considered here we are mostly interested in the case 
where m is small. Unless an integral diverges linearly as m — > one can then drop such 
numerator terms. This is the case for all integrals considered here. 

Let us now discuss how such integrals appear in scattering amplitudes in M = 4 SYM. 
At one loop, the MHV amplitudes are usually represented (or more precisely, their parity- 
even part) as a sum over so-called two-mass easy box integrals [31]. In [13], an alternative 
form was given that uses pentagon integrals with certain twistor numerators. The formula 
given in [13J is a sum over the integrals shown in Fig. |2j Their integrand is given by 
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Xk+1 



"box 
n;i,j 

Figure 2: Integrals contributing to the one-loop MHV amplitude. The box integrand is obtained 
from the pentagon integrand by taking the soft limit p n — ► 0. The dashed line represents here the 
numerator which is denoted by wiggly line in |13j . 



rpent 

n;i,j,k 



pent ji-lijk)(AB(J-ljj+l)n(k-lkk+l)) 



(2.12) 



n;i,j,k (ABi-li)(ABj-lj)(ABjj+l)(ABk-lk)(ABkk+l)' 
We use the notation I for the integrand and F for the function obtained after integration^, 



i.e. 



P n;i,j,k - I 1 n;i,j,k > 



where the number of points n is implicit and appears only in the condition n + 1 = 1. 



A comment is in order here regarding the numerator factors: In [T3] both one-loop and 
two-loop amplitudes are written using numerators with wiggly lines representing 

(AB(i-lii+l)n(j-ljj+l)) := (Ai-lii+l)(B j-lj j+1) - (Bi-lii+l)(Aj-lj . 

. , , , (2 ' 14) 

The MHV amplitude is obtained by changing wiggly lines to dashed lines that stand for 

(ABij), etc. In our discussion we do not distinguish between an integral and the par- 
ity conjugate integral because their difference is a parity odd integral which integrates to 
0(m 2 ). Therefore, in any diagram we can change all (but not just some) wiggly lines to 
dashed lines and vice-versa. In the following we will not distinguish between integrals with 
all dashed lines and those with all wiggly lines. 

We now argue that all integrals appearing in the one-loop MHV amplitudes can be 
thought of as deriving from the pentagon integral. It is clear that pentagon integrals with 
lower number of legs can be obtained from a generic pentagon integral by taking soft limits. 
Therefore we only have to show how to obtain the box integrals of equation (12. 9p . As we will 
see presently, they are also obtained by taking soft limits. Indeed, consider the pentagon 



1 As was mentioned before, some of the integrals are in fact infrared divergent and require a regulator. 
In this section we mostly discuss properties of the (unregulated) loop integrand. It was argued in |13j that 
a correct expression for the loop integrand on the Coulomb branch of N — 4 SYM [7j, up to 0(m 2 ) terms, 
can be obtained by adding masses to propagators on the perimeter of the diagams. 
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jdouble— pent jpenta— box jdoublc— box 

n;i,j,k,£ n;i,j,k n;i,k 

Figure 3: Integrals contributing to the two-loop MHV amplitude. The integrands of the double 
box and penta-box integrals are obtained from the double pentagon integrals by taking the soft 
limits. The dashed line have again the same meaning as in 1-loop case. 



integral of Fig. [2] with k = j + 2, so that only one external leg enters its rightmost corner. 
The soft limit of pj + \ — > corresponds to letting 

Zj+x-xxZj + pZj+2. (2.15) 

Applying this limit to the integrand of the pentagon integral given in (12.121) . we find 

/:;:;::;,., >cw ( 2 - 16 ) 

one reproduces exactly the integrand of the "two-mass hard" box. In this sense we can say 
that the one-loop MHV amplitude is built from a single master integral. We will show in 
section [3] that the latter satisfies a differential equation. 



The situation is almost identical at the two-loop level. The expression given in [13] 
involves only the integrals shown in Fig [3j Recall that we denote the integrand of loop 
integrals by I and the functions obtained after integration by F. For example, the first 
integral shown in Fig. [3] has the following definition, 

j-double— pent 

N (CD(i-lii+l)n(j-ljj+l)) 
nAj "' M (ABCD) (CDi-li)(CDii+l)(CD j-1 j)(CD 

(AB(k-lkk+l)n(e-l ££+!)) 
{ABk-lh){ABkk+l){ABl-ll){ABll+l)' K ' ' 

with the normalisation N — (ijki) and 

pdoublc-pent _ f ^ Z AB d^ZgD ^double-pent /<-, in\ 

r n;i,j,k,t -J in 2 in 2 2 n;i,j,k,e ' 

where the number of points n is implicit and appears only in the condition n + 1 = 1. 
The double pentagon integrals defined in (" 12 . 18[) are in fact infrared finite [13] . Some of the 
pentagon and all of the box integrals we consider are infrared divergent and their definition 
is understood with the mass regulator in place, which leads to the modifications discussed 
above. From the one-loop case it is clear that the penta-box and double box integrals shown 
in Fig. |3]can be obtained from the double pentagon integral by taking subsequent soft limits. 
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x k 



Xl 



Figure 4: Ladder integrals defined in equations (13.190 and (I3.29p . 



So, in summary, at each loop level the integrals appearing in the MHV amplitudes can 
be thought of as deriving from a single master topology. We will show in section [3] that the 
latter integral satisfies a second-order differential equation that reduces its loop order by 
one. 

The relationship between the master integrals and the reduced integrals works at the 
level of the integrand. When the integration is taken into account, in some cases infrared 
divergences can lead to a non-commutativity of the soft limit and the regulator limit. We 
adopt the point of view that in those cases the relation to the master integral still implies 
a certain simplicity of the reduced integral. We show this explicitly for IR-divergent penta- 
box integrals in section [31 which satisfy the same type of differential equations. 

Let us briefly discuss how many functionally different integrals of this class exist. The 
most general double pentagon integral of the type shown in Fig. [3] is the one where four 
of its external legs are doubled (when more legs are added, the function does not change). 
This is possible for the first time at n — 12, e.g. F^^ti^w^ ■ m general, we can consider all 
integrals F^^J pent with zero, one, or two momenta flowing into these corners. Because 
of the symmetries of the pentagon, there are 18 different possibilities. However, because of 
the finiteness of the double pentagon integrals, all of them can be obtained from the most 
general case F^^tfio^ by taking soft limits. Finally, there are 12 penta-box integrals 
F n;Sr b ° X and 6 double box integrals of the type F n d ™ ble " box . The latter are all infrared 
divergent and therefore depend on the regulator m 2 . 



3 Differential equations for loop integrals 

3.1 Finite integrals 
3.1.1 Ladders 

The ladders (or scalar boxes) provide a first example of a class of integrals which satisfy the 
type of differential equations we are interested in. To understand the differential equations 
for the ladder integrals it is convenient to use the dual coordinate notation. For the integrals 
appearing in the two-loop amplitude discussed in section [2] we will pass to the momentum 
twistor notation. We will begin with the one-loop box as an example, using dual variables 
to express the momenta, 
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F^(xi,Xj,x k ,xi) - [ 2 2 2 2 • (3.19) 

J 1,11 -^ir^jr^kr^lr 

Note that the points Xi, Xj, Xk, xi are generic points without any light-like separations. The 
L-loop version of this diagram is shown in Fig. HI Since the integral ( I3.19P is covariant 
under conformal transformations of the x coordinates, it can be expressed in terms of a 
function of the two available cross-ratios, 

2 2 2 2 

•^ij^kl x il x jk / n n\ 

u — ■ ' , V = on - (3.20) 

Thus we have ^ 

-F 1 '- \x{, Xj, Xk, Xi) = ^ ^ , (3.21) 

where the function is known [3"2 1 I3"0" | I3"3"]. In fact the function is best expressed in 

terms of the variables z and z defined byo 

zz 1 

U = — r , V = — r . (3.22) 

(l-s)(l-z)' (l-s)(l-z) V } 

Explicitly it is given by 

$(!)(«,„) = fil ^*f , f (1 \z,z) = log(zz)(Li 1 (z)-Li 1 (^))-2(Li 2 (^)-Li 2 (z)) . (3.23) 

The integral satisfies a simple second order differential equation. The reason for this 
is that acting with the Laplace operator on one of the external points, Xj say, produces a 
delta function under the integral [16], 



D~ = -4ivr 2 5 (4) (x 4 - x T ) . (3.24) 

xj 



This has the effect of localising the integral completely, giving a simple second-order equa- 
tion, 

UiF^^x^x^x^xi) = 2 2 2 - (3.25) 
On the other hand [16], acting on the form of F^ given in ( I3.2ip one obtains 

niPWfax^XtoXi) = ^f^K,M 1] (u,v) (3.26) 

X ik X jl 

where A u „ is a second-order differential operator, 

A UjV = ud 2 u + vd 2 v + {u + v- l)d u d v + 2<9 U + 2d v . (3.27) 



2 Note that z and z are real and independent in Minkowski signature whereas they are complex conjugate 
to each other in Euclidean signature. 
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The equality of the two expressions (13.251) and (13.261) is a second order differential equation 
for the function $W or equivalently Expressing the equation in terms of z and z it 
reads, 

zd z zd 2 fU(z,z) = _£___£_. (3.28) 

The main points we wish to stress are that the action of the operator removes the loop 
integration, leaving a rational function behind. The existence of a simple equation also 
means that the underlying function is a relatively simple pure transcendental function of 
degree two. 

The one-loop integral we have been discussing is just the first in an infinite sequence 
of ladder integrals. These integrals (along with a large class of equivalent integrals [T6] ) 
exhibit an iterative structure, 

/d 4 x r x 2 ji ~r L 
in x ir Xj r x lr 

As before the integrals are invariant under conformal transformations of the Xi and so are 
given in terms of functions of the the two cross-ratios, 



r yXij Xj j X[ ) 2 2" 



(3.30) 



As before the function $ (n ) is best expressed in terms of the variables z and 



& L \u,v)= f(L){ y\ (3.31) 



where 

fV(z,z) = J2 { ~lll L [ 2 r)m )l ^{zz){U 2L . r {z)-U 2L ^.{z)) . (3.32) 

r=0 >' 

The fact that the one-loop box satisfies a differential equation guarantees that all the 
ladder integrals do. One can see that Laplace operator is effectively acting only on a 
one-loop box subintegral which we have already seen reduces to a rational function. Thus 
the operation reduces the loop order of the ladder integral, 

x 2 

D i F ( - L \x i ,x j ,x k ,xi) = -4^^ T F (L_1) (xi, xj, x k , x t ) . (3.33) 

The functions appearing in the explicit expressions exhibit corresponding differential equa- 
tions 

zd x zd s fM{z,z) = fl L - x \z,z). (3.34) 



3.1.2 Pentaladders 

The integrals involving twistor numerators can also satisfy differential equations. We will 
begin with integrals which are similar to the ladders of the previous section. We would like 
to consider the finite pentaladder integrals, beginning with one of the one-loop pentagon 
integrals discussed in section |2j We will write in dual coordinate notation to begin with 
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Figure 5: Figure (a) represents the seven-point integrals defined in (|3.35p and (|3.42p . and (b) 
shows the eight-point integrals defined in (|3.46p and (13.501) . 



and later move to the momentum twistor notation. The integral we would like to consider 
is the seven-point one-loop pentagon integral, which is shown for arbitrary number of loops 
in Fig. |3[a), 

F§\ Xl ,X 3 ,Xt,X 5 ,X 6 )= / . t 2 2 X J T 2 T 2 ■ ( 3 - 35 ) 

Note that we have not normalised the integral so (13.35!) * s dimensionful. The point x a 
is the magic complex point which is null-separated from every point on the null lines 

(x 3 ,x 4 ), (x 4 ,x 5 ) and (x 5 ,x 6 ). It is one of the solutions to 

X a3 = X a4: = X a5 = X a6 = ® ■ (3.36) 

The other solution is its parity conjugate. 

Just like the ladder integrals, the pentagon we are considering is conformally covariant 
so we can write it in the following way 

t pX (Xl, X 3 , X 4 , X 5 , X 6 ) - 2 ~2 ~2 7j " Ia ' ^- d7 ^ 

x 14 x 15 x 36 I 1 " U J 

where u and v are the two non-vanishing conformal cross-ratios, 

2 2 2 2 

u = ^f, v = ^ff. (3.38) 
■^u^e x i5 x 36 

We have chosen to make a factor of 1/(1 — u — v) explicit on the RHS of (I3.37P for later 
convenience. 

In terms of the quantities introduced in section [5] we have 

xlA^F^ ( x i, x 3 , ar 4 , x 5 , x 6 ) = ^ w (u, v) = Ff.^ 5 (3.39) 

The integral satisfies a differential equation very similar to the equation for the ladder 
integrals, 

□i^ = ~ \2 f a x 2 X 2 ■ (3-40) 
This equation translates into a second order p.d.e. for the function 

uvd u d,M l) (u,v) = 1. (3.41) 
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Just as for the ladders we can define multi-loop pentaladder integrals via an iterative 
structure. For example we can consider the integrals, see Fig. EJa), 

F^\x 1 ,x 3 ,x i ,x 5 ,x 6 ) = f ^rf 2 ~T r ^^JnZS)^,^^^- (3-42) 

Equivalently we can write them using the ladder integrals, 

F^\x 1 ,x 3 ,x A ,x 5 ,x 6 ) = [ d ' X 2 r 2 X \ X f 2 F^ L - 1 \x l ,x 3 ,x r ,x 6 ). (3.43) 

As before, we will use conformal symmetry in the x variables to write the integral in 
the form 

F^\ Xl ,x 3 ,x 4 ,x 5 ,x 6 ) = /\ a 2 * (L) ( U > V ) (3.44) 

x 14 x 15 x 36 [L U V) 

Now applying the Laplace operator gives the following differential equation for ty( L \ 

(1 - u - v)uvd u d v ¥ L) (u, v) = (u, v) , (3.45) 

where we define v) = 1 — u — v. 

We can also consider pentaladder integrals with a massive corner on the pentagon subin- 
tegral, see Fig. E(b). At one loop the integral takes the form 

F^( Xl ,x 3 ,x^x 6 ,x 7 ) = f dXr X J . (3.46) 

Here the point x a is light-like separated from x 3 ,X4,Xq,x 7 . It is convenient to write this 
integral as 

\ Via ^ {1) {U,V,W) 

Fy{x 1 ,x 3 ,x 4 ,x 6 ,x 7 )= 2 2 2 ■ r, (3.47) 

p xf^xf^X^j (1 — M — f + UVW) 

where the three cross-ratios are 

22 22 22 

M = V|Z, t, = ^^, u; = ^£. (3.48) 

X U X 37 x 16 x 37 x 36" ( '47 

We remark that the function ^f^(u,v) discussed before can be obtained from ^f^'{u,v,w) 
by taking the (smooth) soft limit w — > 0. 

Following the same logic as before we obtain the equation 

ud u vd v ¥ 1] (u,v,w) = l. (3.49) 

A similar analysis holds for the L-loop case defined by 

/d^x x^ x^ 

. I 2 T 2 7 2 F < ^ L ~ 1 \x 1 ,x 3 ,x r ,x 7 ) . (3.50) 

Writing the integral as a function of the cross-ratios, 

F\'{x 1) x 3) x il x Gl x 7 ) = ■ -, (3.51) 

xf 4 a;f 6 X3 7 (1 — u — v + uvw) 

we find the differential equation, 

(l-u-v + uvw)uvd u d v ¥ L) (u, v,w) = ¥ L ~^ (u, v,w), (3.52) 
with iff®) (u, v, w) = 1 — u — v + uvw. 
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3.1.3 Pentaladders reloaded 



In the case of the pentaladder integrals we can also arrive at the differential equation by 
considering the momentum twistor representation of the integrals. A clue to constructing 
the right operator comes from the fact that the Laplacian naively annihilates the integrand, 

□i— = (naive) . (3.53) 

Of course we have seen that this does not imply that the integral itself is annihilated by 
the Laplace operator because there is an anomaly in the form of the delta function as in 

Let us recall that if one writes the one-loop finite pentagon integral using momentum 
twistors it takes the form, 

*«(«, v) = F 7 Xi 5 = J , (3.54) 

^nt (456[7)(1]234)(AB35) 

7;i,3,5 (AB71)(AB23)(AB34)(AB45)(AB56) ' K ' 

It is possible to construct an operator acting on the momentum twistor variables which also 
annihilates the integrand If.] 1 ^. Let us first introduce some notation to deal with twistor 
derivatives. 

We define a twistor derivative Oy as 

0lJ = Zl 'W^ (3 - 56) 

It acts trivially on four-brackets, Oij(j£km) = iiikm). The normalized integrals are ho- 
mogenous in all external twistors, which can be also written as 

O«J = 0, z = l,...n. (3.57) 

Let us do a trivial exercise which will be important in the following discussion. Let us start 
with the rational function of four-bracket, I = (AB24)/(AB23)(AB34) and act on it with 
the operator 0\%. We immediately get 

{ABU) {AB2A){ABVS) (AB12) 

12 (AB23)(A£34) (AB23) 2 (ABM) (AB23) 2 1 ' ' 

where we used Shouten identity (AB14)(AB23) - (AB2A)(ABl?>) = (AB 12) (AB34). Now, 
we see that the dependence of twistor Z 3 is just through the four-bracket (AB23). If we 
now act with the operator O23, then we get zero, therefore 

°^ pSk) =0 - (3 ' 59) 

Most of the equations we will derive later are based on the same principle. 
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In the case of the pentagon integral the analysis above leads us to define the following 
operator 

d 23 4 = N pl O m 23 N^, jVpi = (456[7)(l]234). (3.60) 
The operator O234 naively annihilates the integrand 1^135, 

0234^7;i°3,5 = (naive) . (3.61) 

We have written 'naive' here because, just as for the Laplace operator, we must remember 
that we are going to perform an integration and there may be anomalies like the delta 
function left behind. In fact we can easily convince ourselves that there are indeed anoma- 
lous terms in place of the naive zero on the RHS of (I3.6ip . Let us express the action of 
the operator on the integral function v). We recall that the cross-ratios take the 

following form in terms of momentum twistors, 

(7123) (3456) (7156) (2345) 

(7134) (2356) ' (7145) (2356) ' 1 ' 

Acting with the operator O on the integral function we find 

6 23 4* (1 '(ti,tj) = - ^^— ) U a^ {UiV) . ( 3. 63 ) 

We already know from acting with the Laplace operator that the RHS of (I3.63P is not zero. 
Indeed the second-order derivative of is 1 according to (13.411) . Writing the equation 
without the factors of N p \ on both sides leads us to the following relation, 

n n [ d ' ZAB (jggg) (7135) 

34 23 J in 2 (AB71)(A£23)(A£34)(AB45)(A£56) (7134) (7145) (3456) ' K ' ' 

This equation can then be used whenever we find the one-loop pentagon as a subintegral. In 
particular one can derive the differential equations for the multi-loop pentaladder integrals. 
As an example we discuss here the finite two-loop pentabox integral, 

^ 2) (u } v) = 

d A Z AB d 4 Z CD (AB35)iV p i(2356) 

~~w 2 m 2 ( AB23) ( AB34) ( AB45) ( AB56) (ABCD) (CD56) (CD71) (CD23) 
Applying the operator O234 and using the identity f)3.64p we find 



(3.65) 



n ) (2356) f d'Zcn N pl (CD35) 

234 {U,V> (3456)7 tn 2 {CD71){CD23){CD34){CD45){CD56)' 1 ' 

Using (I3.63P on the LHS and the definition of the one- loop pentagon (13.5411 . ( 13.551) on the 
RHS we arrive at 

(7135)iV pl . uduVdv ^m^ v) = . (3. 67 ) 



(7134) (7145) (2356) 
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By using cyclic identities one finds that the four-brackets on the LHS can be written in 
terms of u and v and we arrive at the equation 

(1 - u-v)ud u vdM 2) (u,v) = ^ {l) (u,v), (3.68) 

exactly as derived from the Laplace operator in ( 13.45P in the case L = 2. We have rederived 
the equation for the two-loop finite pentaladder to simplify the equations but the derivation 
for L loops is essentially identical and leads to (I3.45P for general L. 

In summary, we have derived the second-order equations (13.451) by acting with the 
twistor differential operator O34O23 on the pentagon sub- integral, thanks to equation (I3.64p . 
What this implies is that whenever we have an integral with the pentagon integral as a sub- 
integral, we can use this mechanism and generate a second-order equation by using f)3.64p . 
Importantly, this is also possible in cases where one cannot apply the Laplace operator, as 
e.g. for the double pentagon integrals that we will discuss in section 13.1.41 Before doing 
this, we are going to discuss the generalisation where the pentagon sub-integral has a mas- 
sive corner. 



In the case of the pentaladders with a massive corner one again finds certain operators 
which annihilate the integrand. Let us consider the integrand of the one-loop pentagon 
with three massive legs, 

.pent _ N pm (AB36) 

8;i,3,6 (AB81)(AB23)(AB34)(AB56)(AB67) ' 1 ' ' 

with N pm = (567[8)(1]234). The integrand is annihilated by certain twistor operators, 

iV pm 24 42 iV-^S i6 = , N pm 75 57 N;^l 6 = . (3.70) 

Applying these operators to the explicit function we actually find that they annihilate it, 
i.e. there is no anomaly associated to this operator on Fg.^Q. In fact we will see that 
Fg^g is given explicitly by |34^l 

^S,6 = * (1) («> v >w) = ^gu\ogv + U 2 (l-u) + U 2 (l-v) + U 2 (l-w) 

-Li 2 (l - uw) - Li 2 (l - vw) , (3.71) 

where the cross-ratios take the form 

_ (8123) (3467) ( (8167)(2356) _ (2367) (3456) 
U ~ (8134) (2367) ' V ~ (8156) (2367) ' W ~ (2356) (3467) ' ^' ' 

We remark that the function ^^(u,v) describing the seven-point integral Ff.^^ can be 
obtained from ^f^(u,v,w) by taking a (smooth) soft limit, 

lim § (1) (m,w,w) = ^ (1) (m,w). (3.73) 

w— >0 



3 Here we are assuming that all dual distances are spacelike. One must carefully analytically continue 
to the regions where some of them become timelike. 
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Using (I3.7ip we then find that the following homogeneous equations are satisfied, 



N pm OM 2 N m \F£™l e = , N pm 75 57 N^Fi™l 6 = . (3.74) 

Writing these in terms of the cross-ratios we find, 

[(1 - u) 2 d u ud u + (1 - w) 2 d w wd w - uv(l - it)(l - w)d u d v 
+ (1 -u){l - w)(l + uw)d u d w -v(l - w) 2 d v d w ]^ {1) (u,v,w) = 0, (3.75) 

and the equation obtained from swapping u and v. As before, we can generalise these 
equations to higher-loop integrals. Note, however, that e.g. the first operator given in 
(I3.74p , O24O42, contains a derivative acting on Z 2 . We can only use (13. 74ft on a pentagon 
sub-integral of a higher-loop integral if the latter does not depend on Z2 elsewhere, to avoid 
cross terms. This is the case if there are at least two additional legs attached at the inter- 
section of the pentagon integral with the rest of the diagram (both additional legs must be 
on the same side of the diagram). 

We will now present a further mechanism for generating differential equations that will 
yield an equation that is valid without the above restriction. To begin with, let us apply 
the differential operator O24O75 to the integrand I^^q- A short calculation shows that 

n n rpent »r (AB36) 

024^75^1,3,6 " "pm ( AB &1) (ABM)* (ABWf ' 

i.e. we have cancelled two propagators. Notice that for obtaining this expression no cyclic 
identities were needed. What this means is that the calculation would go through unchanged 
in the presence of a +m 2 regularisation, so that we do not expect contact terms here. Also 
note that the resulting integral is finite thanks to the numerator factor. Since it is also dual 
conformal, it can only have a very restricted variable dependence. In fact, one can easily 
convince oneself that up to a trivial factor, this integral is a dual conformal three-point 
function that can depend on the dual points X\,x<i,XQ only. Since there are no cross-ratios 
at three points, it must be a rational function. Let us compute the latter. Introducing 
Feynman parameters, we find 

, aia 4 «6(3618) 

ai [aia 4 (8134) + ai« 6 (8156) + a 4 a 6 (3456)] 3 ' 1 ' ' 

The remaining projective integrals are easily carried out, with the result being proportional 
to 

< 3618 > (3.78) 



(8134) (8156) (3456) ' 

Since the differential operator commutes with the normalisation of the integrand we find 

2 a0 75 F^ fi oc N pm - J 3681 ?, . (3.79) 

75 8,1,3,6 pm (8134) (8156) (3456) V ; 

One can easily check this equation on the explicit function f)3.7ip . 
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Figure 6: Six-point double pentaladder integral (|3.8U|) . 



The equation for the sub-integral can immediately be applied to higher-loop integrals 
containing it. There is no subtlety since the twistor derivatives only act on Z 4 and Z$, and 
the latter variables are separated from the rest of the integral. 

We can write down several other operators which follow the same pattern. The idea is 
simple: if we can manage to reduce the number of propagators to three, dual conformal 
symmetry will imply that the resulting integral is rational]^ An example is given by O42O75. 
Note, however, that in that case the corresponding equation can only be used iteratively in 
cases where there is no Z2 dependence in the remainder of the integral. 

In conclusion, whenever we encounter a pentagon sub-integral of the type considered 
in this section, either with or without massive corner, we can always find at least one 
differential equation that can be used iteratively. In some cases we can even find further 
equations. In particular what this means is that we can find differential equations for 
all double pentagon integrals appearing in the MHV amplitudes. In the next section we 
are going to spell these equations out explicitly for the case of an infinite class of double 
pentaladder integrals. 



3.1.4 Double pentaladders 

The fact that we have shown that one can derive differential equations using the momentum 
twistors is important because it allows us to find equations in cases where no Laplace 
operator derivation exists. As an example we will consider the six-point finite double 
pentaladder integrals. These integrals can be defined in terms of the pentaladders of the 
previous subsection, 

/d A x x 2 x 2 
2 2 a \ H Fpi~ 1 \x r ,x 4 ,x 5 ,x 6 ,x 1 ) , (3.80) 

see Fig. |6j Conformal symmetry in the x variables tells us that we can write 

x 2 

F^'(x 1 ,x 2 ,x 3 ,X4,X5,x 6 ) = 2 f 2 Vt {L) (u,v,w) , (3.81) 

X 26 X 35 X 14 

4 Naively one might find it surprising that dual conformal symmetry is helping us here, since the momen- 
tum twistors already made that symmetry manifest. Indeed with six twistor variables it is possible to form 
homogeneous cross-ratios but the final result contains no extra functional dependence on these variables. 
The reason is that we are combining the symmetry with the fact that our functions are Feynman integrals 
built from propagators, which leads to additional restrictions on the possible dependence on the twistor 
variables. 



19 



where u, v and w are the three cross-ratios, 



2 2 2 2 2 2 

M = ^|6, v = ^ff, w = ^ff- (3-82) 

X 14 J '36 x 14 x 25 x 25 x 36 



To derive differential equations for these integrals we need to express them in momentum 
twistor notation. We will consider the finite six-point two-loop double pentagon integral as 
an example, 

rt 2 \u,v,w) = F^ 4 y™ t (3.83) 
d 4 Z AB d 4 Z CD {AB13){CDAQ)N dp 



VK 



2 m 2 (AB6l)(ABl2)(AB23)(ABM)(ABCD)(CD45)(CD56)(CD6l) ' 



where N dp = (2345) (5612) (3461). In terms of the momentum twistors the three cross-ratios 
take the form, 

_ (6123) (3456) _ (1234) (4561) _ (2345) (1256) 
U ~ (6134) (2356) ' V ~ (1245) (6134) ' W ~ (2356) (1245) ' ( } 

We see that the double pentagon integral contains the one-loop pentagon that we studied 
in the previous subsection as a subintegral. Thus we can apply the same kind of differential 
operator and derive a differential equation. The operator we wish to apply here is 

612 = N^OuO^N^ 1 . (3.85) 

Acting with this operator on the the double pentagon integral and using the relation (I3.64p 
with the appropriate relabelling we find 

1 r d 4 Z CD N dp (CD13)(CD46) 
~ "(1234) J iir 2 (CD12) (OD23) (OD34) {CD Ah) (CD5Q) {CD61) 

(3461) fi«( M ,^). (3.86) 



(1234) 

Here Q^'(u,v,w) is the one-loop finite hexagon integral, 

K ' ' ' J m 2 (CD12)(CD23)(CD34)(CD45)(CD56)(CD61) ' V ' ; 
where N hex = (2345) (5612). 

From the discussion in section 13.1.31 it is clear that the above equations can be straight- 
forwardly generalised to the case where the pentagon sub-integrals have a massive middle 
leg. 

3.2 Divergent integrals 

Some of the integrals appearing in the MHV amplitudes are infrared divergent. These 
divergences are regulated by masses as explained in section [2j In this section we show that 
one can find similar differential equations as above, that hold up to 0(m 2 ). 
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3.2.1 One-loop box integrals 

As explained in section 121 they are special cases of the general pentagon topology with 
magic numerator, when no external legs are attached to the rightmost corner. In this case, 
one obtains a two-mass hard box integral, whose integrand is given by 



,-i.ox I d A Xp 44 |"» ss> 

n ' 1 ' 2 I m 2 (x 2 + m 2 )(x 2 + m 2 )(x 2 + m 2 )(x 2 + m 2 ^ ' 



Unlike the generic pentagon integral F^- k the two-mass hard integral is infrared divergent 
and depend explicitly on the mass m that regulates the infrared divergences, see section 
|2j Here we discuss that the differential equations we found for the finite integrals work in 
exactly the same way here, up to 0(m 2 ) corrections. To understand why this is the case, 
consider again the Laplace operator, but acting on a massive propagator^! We have 

= ~ 8 ( 2 ""l 2 , 3 =-4nr 2 ^)K t ) + (9(m 2 ), (3.89) 
4 + m l (4 + m 1 ) 6 

so that in this case we would indeed "undo" one integration as in the massless case, but up to 
0(m 2 ) corrections. As we discuss presently, this leads second-order differential equations for 
this integral, at least for n > 5. We will first discuss the generic case where ( 13.881) represents 
indeed a two-mass hard integral, and subsequently discuss the degenerate one-mass case 
where i = 4 or i = n. 

In the former case, there is one off-shell leg that we can apply the Laplace operator to, 
see equation (I3.89[) . This leads to a differential equation valid to 0(m 2 ), namely 

444 □* (44) -1 *5Ta = ~ 4 + °( m2 ) • ( 3 - 9 °) 

We can explicitly verify (I3.90p by differentiating the known answer for the two-mass hard 
integral. In terms of Feynman parameters, it is given by 



| 4 4 $( a 2 + a 3 + a 4 + Oii - 1 



F b ° X = / dOii X 2^3i "V^2 TU3Tu jT u,- v 

n ' r ' 2 J 1 \oliol±x\± + ai(a 2 xl im + a 3 x| im + «44,m) + m2 ] 2 ' 
where 4 m '■= 4 + m2 > e ^ c - One obtains 

F^2 = \ log 2 (2/12/22/a) + 2 Li 2 (1 - y 2 ) + 2 Li 2 (1 - y 3 ) + 0(m 2 ) , (3.92) 

where we used the variables yi = x^Jrv? ,y 2 = 4/ a; 2i ;2/3 = x \%l x \%- J us ^ hke in section 
13.11 we can also find twistor derivatives that lead to differential equations. However, the 
two- mass hard integral is infrared divergent, and here we have dropped certain 0(m 2 ) 
terms. When applying differential operators to an infrared divergent subintegral, e.g. in 
the double box case, one must be very careful when making such approximations. We 
postpone a detailed analysis of the double box integrals to future work. 
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Figure 7: Five-point penta-ladder integrals considered in the text. 



3.2.2 Penta-box integrals 

A class of finite pentaladder integrals was already discussed in section 13.11 Here we show 
that the IR divergent penta-box integrals also satisfy analogous differential equations, up 
to 0(m 2 ). 

Let us consider the pentaladder integrals shown in Fig. [7J The discussion for other 
divergent pentaladder integrals would be very similar. At one loop, we have 

t f d 4 Z AB (2345) (4512) (^13) 

5;5 ' 4 ' 3 J m 2 {51AB)(12AB)(23AB)(MAB)(45AB) ' 1 ' 

Since all external legs are on-shell, we cannot use the Laplace operator. We already know 
from section 13.1.31 that in such cases twistor derivatives can have the same effect. Indeed, 
consider the action of the differential operator 

"^^^■■m^-^&'-k- (3 ' 94) 

where N = (2345) (4512), One can check that it annihilates its integrand, up to 0(m 2 ) 
terms. As discussed above, one can think of the mass as regulating distributional terms 
that arise when acting with the derivatives. Let us now find out what happens when we 
act with the operator in (I3.94p on the function obtained after integration. Computing the 
one-loop integral F^i 3 directly using Feynman parameters in the small m 2 limit, we have 
(see Appendix [B]) 

Fg^ts = -^log 2 \yiym) - 2Li 2 (l - y x ) - 2Li 2 (l - y 2 ) + ^ + 0(m 2 ) . (3.95) 

Here y\ = x 35 /xf 3 ,y2 = x\ h jx\^ y 3 = x\Jm 2 y and we recall the relation of the dual variables 
to the momentum twistors, 

_ (5123) _ (1234) _ (2345) _ (3451) _ (4512) 
13 (51)(23)' 24 (12)(34)' 35 (23) (45)' 14 (34) (51) ' 25 (45) (12) ' 1 } 



5 In this case one could also note that massive propagators are formally identical to free AdS bulk-to- 
boundary propagators, and use the corresponding equation of motion, see [35j . For this it is necessary to 
consider situation with generic masses rrn 1 but this is no restriction thanks to dual conformal symmetry. 
This would yield equations that are exact in the rrii, i.e. that do not rely on taking the masses small w.r.t. 
the Mandelstam invariants. It is an interesting open question how to extend this efficiently to the on-shell 
case. 
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Using the chain rule and acting on (13.95!) one can easily verify that 
(5123) 



(4513) 



d „ d , (5134) d 



= "I + 0(m 2 ) , (3.97) 



<9Z 2 <9Z 3 (1245) dZ 3 _ 

i.e. the r.h.s. is not zero but a simple rational term. In other words, it is still true that the 
differential operator "undoes" the integration, at least up to 0{m 2 ). 



We now discuss the all-loop generalisation of these equations. The two-loop version of 



IS 



^penta-box f d*Z AB d 4 Z CD (j7T 2 )~ 2 N(CD13) 

where N = (3451) (2345) (4512). 

As in the previous sections, we need to discuss the pentagon subintegral. In fact the 
subintegral is the same finite integral that was discussed in section 13.11 with the difference 
that since the two-loop integral is IR divergent its propagators contain +m 2 pieces. How- 
ever, we already know that the special numerator makes the subintegral finite, and therefore 
the effect of the masses in the pentagon propagator should be only 0(m 2 ). Therefore 
we can neglect these masses when acting with the differential operator, and it is clear that 
we will obtain differential equations completely analogous to those discussed in section 13.11 
The only difference is that by construction they are only valid up to 0(m 2 ) corrections. 
Therefore we have 

jUjiVO^iV- 1 F:V^ - = Jg*. + 0{m 2 ) . (3.99) 

Of course one can immediately generalise this equation to arbitrary loop order, by replacing 
the integral on the LHS and RHS by the L-loop and [L — l)-loop penta-box integrals, re- 
spectively. Finally, the other divergent penta-box integrals which differ in the arrangement 
of external legs all satisfy analogous differential equations. 

In summary, we have found differential equations for all integrals contributing to the one- 
and two-loop MHV amplitudes in M = 4 SYM, except for the six double box integrals. The 
differential equations we find are iterative in the loop order and apply to infinite series of on- 
shell loop integrals. We expect that integrals contributing to higher-loop MHV amplitudes 
and other helicity configurations are constrained in a similar way. 

The same mechanisms for generating differential equations that we described above 
should also hold for the double box integrals. However a separate analysis is required there 
since they do not contain IR finite subintegrals. We leave this for future work. (An en- 
couraging fact is that in Appendix |A] we find differential operators that naively annihilate 
the double box integrand.) Of course, the double box integrals are in some sense the sim- 
plest integrals appearing in the two-loop MHV amplitudes, so perhaps they could also be 
computed by other means. Also, they are related to the more general integrals that do sat- 
isfy differential equations through soft limits. If one can control the non-commutativity of 
the soft limit with the regulator limit this can also provide a useful tool for computing them. 



6 This argument presupposes that the divergent box subintegral produces only logarithmic but no linear 
divergences in to 2 . This is the case for all integrals we consider. 
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4 Solving the differential equations 



4.1 Harmonic Poly logarithms 

We will now discuss solving the differential equations. The solutions we find are typically 
written in terms of polylogarimic functions. Here we will give a brief introduction to the 
notation we use which can be found in [361 [37] . 

The simple harmonic polylogarithm functions are given in terms of iterated integrals. 
We start from the functions 

#(0;x) = logz, #(l;:r) = -log(l-x). (4.100) 

Their derivatives are given by the fractions 

^H(a;x) = f(a;x) } f(0;x) = -, /(l ; x) = _L_ . (4.101) 
ax x 1 — x 

We can now define weight w harmonic polylogarithms depending on a vector with entries 
of ones and zeros. For a vector of all zeros of length w (written W ) we have 

H(0 w ;x) = ^\og w x. (4.102) 

Further harmonic polylogarithms depending on a more general vector m w = (a,m ffi -i) can 
be defined by repeated integration, 

H(m w ;x)= / dx'f(a; x')H(m w _ 1 ; x') . (4.103) 
Jo 

Below we save space by writing the vector subscript. Finally we can use the 

shorthand notation where proceeding from right to left a zero is eliminated if it is to the 
left of a non-zero entry while at the same time one is added to the value of that entry. For 
example 

#3,2(a0 =#(0,0, 1,0,1; re). (4.104) 

There also exist harmonic polylogarithms with negative arguments but we will not need 
them here. 

The simple harmonic polylogarithms defined above are all single variable functions. 
However they can be generalised to functions of more than one variable. In general the 
relevant functions are the Goncharov polylogarithms. Here we will only need functions 
of two variables at most which can be described in terms of two-dimensional harmonic 
polylogarithms. These functions are defined in much the same way as the simple harmonic 
polylogarithms but we allow slightly more general arguments of the vector m w . This is 
done by enlarging the set of weight one functions and fractions to include 

#(1 -y;x) = - logfl - -1—) , /(l -y;x) = - . (4.105) 

V 1 — ys 1 — x — y 
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4.2 Solution for seven-point pentagon and penta-box integral 

Let us consider the equation for the one-loop pentagon integral. The equation is 

uvd u d v ¥ 1 \u,v) = l. (4.106) 

The general solution is 

^ {1 \u,v) = log u log v + g(u) + g(v) (4.107) 

for some single variable function g. 

Let us recall that we chose to write the pentaladder integrals is the following way, 

fJ 1 L) (xi,X3,X4,x 5 ,x 6 ) = %la — ^' V \ . (4.108) 

This form makes it clear that the function ^(^(m, t>) should have a simple zero when 
1 — u — v — since the integral clearly has no special singular behaviour on this region. 
Thus we have 

¥ L \u,l -u) = 0. (4.109) 
The boundary condition (14. 109j) implies that the function g should satisfy 

g(u) + g(l -u) = -logwlog(l - u) . (4.110) 

The solution of (I4.110p is ambiguous up to any function odd in swapping u and 1— u. Thus 
we find 

g(u) = Li 2 (l -u)- |C(2) + h(u) - h(l - u) , (4.111) 

for any single variable function h. 

We need some extra data to fix the function h(u) — h(l — u). If we adopt maximal 
transcendentality (which translates here to the fact that h(u) can only be a sum of harmonic 
polylogarithms) and also the requirement that there are no branch points of ^/^(u,v) for 
u > 0, v > then we fix h(u) = 0. Indeed if we assume that we can write h(u) as a sum of 
harmonic polylogarithms of degree 2, 

h{u) = aLi 2 (w) + b log 2 u + c logwLi^w) + dU\{u) , (4.112) 

then the only solution where h(u) — h(l — u) has no branch points at u — 1 is indeed h = 0. 
Thus we deduce that 

^ (1) (u, v) = log u log v + Li 2 (l -u) + Li 2 (l - v) - C(2) . (4.113) 

Given the correct answer for we can go on to determine the next pentaladder function 
\I/( 2 ). The equation we have to solve is 

(1 -u-v)uvd u d v ^ i2 \u,v) = ^ l) (u,v). (4.114) 

It is helpful to divide by the factor (1 — u — v)uv and use partial fractions to write 

d u d v ¥ 2 \u,v) = 



uv 



u(l — u — v) ^ v (1 — u — v) 



(4.115) 
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Now we look for particular solutions which reproduce each of the three terms on the RHS, 
i.e. 

d u d v H°(u,v) = ^-J, (4.116) 

uv 

d u d v H u (u,v) = (4.117) 
u{l — u — V ) 

d u d v H v (u,v) = ^^-. (4.118) 
t>(l — u — V) 

We find the following solution for H°, 

H°(u, v) =2C(4) + i((2) 2 - C(2) (Li 2 (l -u) + Li 2 (l - „)) 

+ logMlogwdlog-ulogt; + Li 2 (l -u) + Li 2 (l - v ) - C(2)) 
+ log«(2# 2il (l -v)- Li 3 (l -u) + C(3)) 
+ log«(2/T 2 ,i(l -u)- Li 3 (l -v) + C(3)) 

- 3(^,1(1 - «) + #3,i(l " «)) , (4.H9) 

where we have used the notation of harmonic polylogarithms [36]. In fact ii/2,1 is a Nielsen 
polylogarithm and can be reexpressed in terms of ordinary trilogarithms but we prefer not 
to do so here. For H u we findll, 

H u (u,v) =C(2)(#o,i-„(«) - Li 2 («) - log«log(l - u)) + H 0fi (u)H lfi (v) 
+ \ogv(H 0A ^ vfi (u) - H 2 fl(u)) + loguFi )liQ (u) 

- #2,1,0(11) + #o,i-«,i,o (w) + H 0>1 _ v (u)Hi fi (v) , (4.120) 

where in addition to harmonic polylogarithms we have used the two-dimensional harmonic 
polylogarithms of [3T|. We have neglected some obvious simplifications such as H fi(u) = 
|log 2 w and Hi t0 (u) = — Li 2 (u) — log it log (1 — u). By symmetry we have also 

H v (u,v) = H u {v,u) . (4.121) 

Thus the general solution to the differential equation (14.1141) is 

^ (2) (u, v) = H°(u } v) + H u (u, v) + H v (u, v) + g{u) + g(v) , (4.122) 

for some single variable function g. The function g can be constrained by the fact that 
iffW(u, v) should vanish when 1 —u — v = 0. Again assuming that it is purely transcendental 
and has no spurious branch cuts (although as presented, the branch cut structure of H u is 
not obvious) fixes the function. We find 

g{u) = - K(2) 2 - 3C(4) + C(2) (| log 2 u + Li 2 (l - «)) (4.123) 
+ logu(Li 3 (l -u) + # 2 ,i(l - u) - 3C(3)) + 2fT 3jl (l - u) + 3# 2>lil (l - u) . 

We have checked that this combination of functions numerically matches the Mellin-Barnes 
representation of the integral in the triangular region u > 0, v > 0, u+v < 1 to a high degree 
of accuracy. Thus we are confident that (14.1221) together with fl4.119p . fl4.120p .f l4.12ip . and 
(I4.123P correctly give the two-loop pentaladder function. 



We thank L. Dixon for pointing out a typo in a previous version of this formula. 



26 



4.3 Equations for eight-point pentagon and pentaladder integrals 

We would now like to discuss the equations for integrals with pentagon subintegrals with 
massive corners. For example we have the pentaladder integrals with a massive corner, 

p( L ) ( T T T T t n x la ¥ L \u,v,w) 

r p i {Xi,X3 J X4,x G ,x 7 ) — — — 2 — 2~7i : r> (4.124 J 

1 xf 4 xi 6 x^ 7 (1 — u — V + uvw) 

Note that the way we have chosen to write the RHS means that ^( L ) (u,v,w) vanishes when 
1 — u — v + uvw = because the integral has no pole in this region. 

As we have seen for the pentaladders with a massive corner we have two sources of 
differential equations. Firstly from applying the Laplace operator we have: 

(1 - u -v + uvw)uvd u d v ¥ L) (u,v J w) = ¥ L - 1] (u, v, w) . (4.125) 

Then in the twistor language we can apply the operator 

24 75 (4.126) 

which also produces a rational function. In fact we have 

N vm 2 4P 75 Nll¥ 1 Hu,v,w) = 24 75 ¥ 1 \u,v,w) = - N ^ 368 y . (4.127) 

pm 75 pm v > > ) 24 75 \ , , j (3456) (3481) (5681) V ; 

This leads to the equation 

w[(l - w)d w ((l - w)d w + v(l - v)d v + u(l - u)d u ) + u(l - u)v(l - v )d u d v ]^/ {1 \u, v , w) 
= 1 - u - v + uvw = ^ (0) (u,f,w) . (4.128) 

More generally we have for the L-loop pentaladder with a massive corner, 

w[(l - w)d w ((l - w)d w + v(l - v)d v + u{l - u)d u ) + u{l - u)v(l - v)d u d v ]¥ L \u, v, w) 
= ¥ L - 1 \u,v,w) . (4.129) 

From (I4.125P we know the result of acting with the final term in the operator on the LHS 
so we can rewrite the equation in a simpler way. We find 

wd w [(l - w)d w + w(l - u)d u + v(l - v)d v ]¥ 1) (u,v, w) = 1 - u - v , (4.130) 

or more generally for L loops, 

wd w [(l-w)d w +u{l-u)d u +v(l-v)d v ]¥ L \u,v,w) = (i-^-^)^ ^MM^ _ (4131) 

1 — u — V + uvw 

Now let us consider (14.1251) restricted to w — 1. We find that the operator factorises, 

u(l -u)d u v(l -v)d v ¥ L) (u,v, 1) = ¥ L - 1) (u,v,l) . (4.132) 
In the one-loop case we then have simply 

uduvd^iu, v, 1) = 1 . (4.133) 
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and so 

^ {1 \u,v,l) = logwlogu + f(u) + f(v) , (4.134) 

for some single variable function /. We recall that ^^ L \u, v , w) vanishes when 1 — u — v + 
uvw = 0. In our case where w = 1 this implies that ^^'(u,v, 1) vanishes when u = 1 or 
v — 1. Thus we have 

/(u) + /(l) = 0, (4.135) 

which implies that /(tt) = and so 

^ (1) (m,v, 1) = logulogv. (4.136) 
In the two-loop case it is again simple, 

u(l-u)d u v(l-v)d v ¥ 2 \u,v,l) = log u log v . (4.137) 

Integrating we find 

¥ 2 \u, v, 1) = [H 0i0 (u) + H 1)0 (u)][H 0fi (v) + H 1>0 (v)) + f(u) + f(v) (4.138) 

for some single variable function /. As before we have the condition that ty( L \u, v , 1) 
vanishes when u = 1 or v = 1. Thus we find 

- C(2)[#o,o(u) + #i,o(«)] + f( u ) + /(!) = • ( 4 - 139 ) 
Setting u — 1 this tells us 

/(1) = -|C(2) 2 . (4.140) 

Thus we have 

f(u) = C(2)[H , {u) + H lfi {u) + |C(2)] (4.141) 

and we have solved for ^^(u^, 1) quite simply. 

Now we can return to the other equation (14.131)) . Let us perform the w integral, 

[(1 - w)d w + u(l - u)d u + v(l - v)d v }^( L) (u, v, w) = X {L) (u, v, w) + f£\u, v) (4.142) 
where 

X^(u,v,w) = (l-u-v) / (4-143) 

J 1 t{l — u — V + uvt) 

and fx{u, v) is an arbitrary two- variable function. We can determine fx however by setting 
set w — 1 in ( 14. 142)) and we find 

Kl -u)d u + v(l -v)d v ]¥ L \u,v, 1) = f ( x L \u,v) . (4.144) 

The LHS is known from the Laplace operator as shown above so given ^^(-u, v, 1) we know 
fx (u,v). In the one-loop case we have 

f£\u,v) = (1 -u)logu + (l -«)logu. (4.145) 

while 

f w dt 

X (1 \u,v,w) = (1 -u-v) / — = (1 -u-v)\ogw. (4.146) 

Thus (14.1421) becomes a first order equation for ^^'(u, v, w). One can easily convince oneself 
that the unique solution is ( 13.711) given the known boundary condition for v, 1) or 

equivalently ^ w (u,v,0) = ^ w (u,v). 

More generally we always have a first order equation for ^( L \u, v, w) and thus one can 
write the function in a simple integral form. 
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4.4 Equations for the six-point hexagon and double-pentagon in- 
tegrals 

Let us now discuss the simplest cases of the double pentaladder integrals. To start we will 
consider the finite six-point double-pentagon, defined in ( I3.83p . In terms of the twistor 
variables the independent cross-ratios take the form 

(6123) (3456) (1234) (4561) (2345) (1256) 

(6134)(2356) ' (1245) (6134) ' (2356)(1245) ' { ' ' 

Translating the equation f l3.86|) into an equation for the function we find 

(l-u-v + ^)uvd u d v + U vwd v d w + ^ V uwd u d w - ^-d w (l - w)d w Q {2) (u,v,w) 
= n {1 \u,v,w). (4.148) 

with e = y/5 2 — Auvw + 5 and 6 — u + v + w— 1. 

In fact the sign of the square root is ambiguous and so the equation separates into two 
second order equations, one homogeneous, the other inhomogeneous, 

(l — u — v) uvd u d v — uvwd v d w — uvwd u d w QS 2 \u, v, w) = 0^(u, v, w) , (4.149) 
uvd u d v + uwd u d w + vwd v d w — wd w (l — w)d w Q^ 2 '(u, v, w) = . (4.150) 

In fact we can set w = in ( 14.149)) and we find 

(1 -u-v)uvd u d v tt i2) {u,v,0) = tt {1 \u,v,0) . (4.151) 

This equation is identical in form to that for the functions ty( L \u,v) from the pentaladder 
with a massless corner. At one loop we have in fact 

fi (1) (w, v, w) = log w log v + Li 2 (l -u) + Li 2 (l - v) + Li 2 (l - w) - 2((2) , (4.152) 

and indeed n' 1 '^,^, 0) = (u,v). We have verified numerically from the Mellin-Barnes 
representation that the same holds at two loops, Q( 2 \u,v,0) = ^^ 2 \u,v), so this two- 
variable slice of the function is the same function that we have solved for from the 
pentaladder equation. 

Combining the two equations (I4.149|4. 150)) we can write a factorised operator, 



-u 



(l-u)d u -v(l-v)d v + (l-u-v)(l-w)d w ]n {2 \u,v,w) = fi (1) (w,t',w) . (4.153) 



After performing one integration on w we have a first-order equation for the function £1^. 



5 Conclusion and outlook 

In this paper we derived a new type of differential equations for on-shell integrals. The 
equations are second order and reduce the loop order by one. We presented several classes 
of infinite of integrals closed under the differential equations. In particular, these include 
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Figure 8: The iterative differential equations can always be obtained if a key sub- integral is 
present, as shown in the Figure. The salient feature is the presence of a "magic numerator", as 
indicated by the dashed line. The number of differential equations that can be obtained for an 
integral depends on the precise configuration of external legs. 

most integrals appearing in M = 4 MHV amplitudes to two loops. One may speculate that 
the same should be true to all loop orders. 

We also discussed how such differential equations can be solved by making reasonable 
assumptions about the boundary conditions. The explicit solution for a finite two-loop 
seven-point and a divergent two-loop five-point integral was given. The answers are rela- 
tively simple and one can hope that one can solve for infinite series of such integrals, as 
is possible for the off-shell ladder integrals [32j [301 E3] . This would be an important step 
toward summing the perturbation series in Af = 4 SYM. See [3S] for a study of the sum- 
mation of the off-shell latter integrals. 

The main motivation of this study was to contribute to the determination of the all-order 
S-matrix of M = 4 SYM. One may ask whether there is a short-cut by deriving differential 
equations for the whole amplitude, as opposed to individual integrals. In Appendix [A] we 
found differential operators that annihilate (naively) the integrand of the one-loop MHV 
amplitude. In order to understand these differential operators better it is important to 
systematically investigate distributional terms that can appear and lead to inhomogeneous 
terms in differential equations. A related open question is the link of our differential equa- 
tions to the conformal symmetry of M = 4 SYM. 

We expect our differential equations to be useful for other theories as well. Due to 
their nice properties it is certainly desirable to include the integrals with magic numerators 
introduced in [13] and discussed in [27] and the present paper in the integral basis. The 
necessary numerator identities were given in [27]. Moreover, one can find differential equa- 
tions of the type presented here for larger classes of integrals. This is clear from the way the 
differential equations were derived: all one needs is a certain sub-topology with "magic" 
numerators. In the cases considered in this paper the latter appeared on top of a pentagon 
topology, see Fig. [BJ which made sure that the integrals are dual conformal. However, one 
could drop this requirement and extend the current method to integrals that are not dual 
conformal. 

Another natural extension of our work is to use different regulators. Some of the in- 
tegrals discussed here are finite, and in the case of infrared divergences we used a mass 
regulator, allowing us to stay in four dimensions, which is very natural when working with 
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twistor variables. However, this is not mandatory. Although the details may be slightly 
different, the mechanism for finding differential equations should work with other regulators 
as well. E.g. in dimensional regularisation one might expect similar equations to hold, at 
least to some order in e, where the dimension is 4 — 2e. 

It is interesting to obtain equations exact in m 2 , which would allow to discuss e.g. 
scattering amplitudes on the Coulomb branch of M = 4 SYM. Equations exact in m 2 
could also be useful to study integrals appearing in realistic amplitudes involving massive 
particles, see e.g. (39]. It should be possible to derive such equations at least in some cases, 
perhaps by noticing that four- dimensional massive propagators are formally identical to 
AdS bulk-to-boundary propagators [35] . 
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A Twistor differential operators 

In section El we showed that all integrals appearing in the one- and two-loop MHV ampli- 
tudes with at least five external legs satisfy differential equations. In order to obtain these 
equations it was crucial to find differential operators that can act on a sub-integral and 
reduce it to a rational function. 

In order to find the relevant twistor differential operators we found it useful to system- 
atically write down all differential operators that have of property of (naively) annihilating 
the integrand of a given integral. Of course this is not equivalent to finding a useful differ- 
ential equation for the function after integration, but it served as a helpful starting point 
for our above analysis. Here we present more details of the mechanism how differential 
operators annihilate the integrand of the integrals discussed above. The reason is two-fold: 
firstly, it would be interesting to understand in more detail how the distributional term 
that freezes the loop integration is generated. Secondly we expect our differential equation 
method to apply to larger classes of integrals. We hope that this section will prove a useful 
starting point for such investigations. 

We do not claim that we exhaust all differential equations for a given diagram, but 
we show that there is always at least one for all diagrams appearing in the one-loop and 
two-loop MHV amplitudes. 

A.l Diagrams in the one-loop amplitude 

The form of one-loop MHV amplitude presented in [13] contains sum of two- mass-hard boxes 
and double pentagons with "magic numerator". As was mentioned above, for convenience 
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we use the pentagons with dashed lines, the difference being a parity odd piece that vanishes 
after integration (or is of order 0(m 2 ) in the massive regularization). 



A. 1.1 Two-mass hard diagram 

The two-mass-hard boxes as integrals in momentum twistor space are^l 

This integrand is annihilated trivially by the first order differential operators O21, O34, 
O u +i, Oi-\- U and also for all k and j — 5, . . . i — 1, i + 2 . . . n. 

We are interested in obtaining non-trivial second order differential equations for /^m/r 
Applying the operator 2 4 on the un-normalised integrand hmh/N we get 

1 1 
° 24 (ABU) (AB23) (ABM) (AB ii+1) = (AB12)(AB34) 2 (AB ' ( A - 155 ) 

and this is trivially annihilated by Oyi and O43. Similarly, we can start with the operator 
O31 and then apply 0\ 2 and 43 . As a result, we get four differential equations^ 

Ol224hmh — O^ihmh = Ol231^2m/i = O^i^mh = , (A. 156) 

where 

O m = NO^OmN- 1 . (A. 157) 

In the boundary case of i = 4 or i = n (which corresponds to a one-mass diagram) we 
cannot use both "directions" of the derivatives. So in the first case we use only O31 to 
create a massive corner, while in the second case only O24. There is a special case of n = 4 
which corresponds to the zero-mass diagram. In this case we cannot use our construction. 
However, this diagram itself represents the one-loop MHV amplitude and we will see at 
the end of this subsection that there exist differential operators that annihilate naively all 
MHV one-loop integrands (including the n = 4 case). 

For reasons that will become clear later, it will be quite useful also to find a differential 
operator acting only on a limited number of twistors. Suppose that the only operators we 
want to use are 0^2, Ok3, O n i> O54 for any k. Even in that case there exists a second-order 
operator, although somewhat complicated, 



~ (2345) - (2456) ^ < (1356)(2345) ^ ' 

° 4354 _ (1234) ° 1223 _ (2356) ° 54 + (1234) (2356) ° 23 



l2mh 



0, (A.158) 



and also its flipped version. It is hard to find an explanation for this operator (like we did 
for other operators above), but its existence is important for the discussion of the two-loop 
double boxes. 



8 We changed the labels from i — > i + 1 w.r.t. to the main text. 

9 Note that we do not include the second order operators that include O21 and O34 as first order pieces. 
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A. 1.2 Pentagon diagram 

The pentagon integrals in momentum twistor space are 

T ij _ N(AB(m)n(ir-iu+i)) 

pent (AB12)(AB23)(ABi-li)(ABii+l)(AB]]+l)' 1 ' 

where N = (2ijj+l). The strategy is very similar to the two-mass hard case. We can 
find many trivial first order operators that annihilate the integrand, and we can use again 
the same trick of creating a massive corner. Here we have even more freedom in the choice 
which propagator we want to kill, e.g. we can first apply the operator O13, yielding 

(AB (123) n (i-lii+1)) 



( )l3 ( AB 1 2) ( AB23) ( AB i-li) ( AB ii +1 ) ( AB j j +1 ) 

(AB(123)n (A16Q) 



(AB23) 2 (AB i-1 i) (AB i i+l) (AB j j +1) ' 

which is annihilated by the operator O32. Finally, we can find four operator^] that non- 
trivially annihilate Ip en ti 

O&izlpent = Oi22,lIpent = ) (A. 161) 

and the same for (1,2,3) — > (i+l, i, i-1). We see that the operators act separately just on 
the left or right region. 

Now, something interesting happens when i = 4, i.e. when the top corner is massless. 
Then both regions (1,2,3) —> (i+l, i, i-1) overlap and the equations we wrote are not 
valid anymore. However, the numerator saves the situation and helps us to find another 
differential equation. First, let us apply the operator O23, 

(AB(123) n (345)) 



O 



2:-! 



(AB12) (AB23) (ABM) ( AB45) (AB j j +1) 

(AB(123) n (245)) (AB24) (AB( 123) n (345)) 



(AB12)(AB23)(AB34)(AB45)(AB jj+1) (AB12)(AB23)(ABM) 2 (AB45)(AB jj+1) 
( 1234 ) 

~ (AB12)(ABM)(AB j j+1) ' ( ' 162) 

where we used (AB(123) n (234)) = (AB23)(1234). The term we obtained is annihilated by 
0\2 and O34. Similarly, we can start with the operator 43 and repeat the above calculation. 
In total we get four differential equations, 

Ovmlpent — du23lp J ent = ^5443^^ = ^3243^pfnt = • (A. 163) 

Here there are no boundary cases for j = 5 or j = n. 

In addition to these operators, we can also find two more special operators that annihi- 
late the integrand of the pentagon. Like in the two-mass hard case they also contain only 
a restricted number of derivatives, but their forms are much nicer, 



i-2 



Oki-nk + On—I 



k=2 



± pent 



Ok3i+lk + O^i+i 



k=A 



4it = 0. (A. 164) 
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For the boundary case when j = i + 1 or j = n some operators are not present. 
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We can see that two of equations for i = 4 were exactly of this type. 

We completed the picture of differential operators for all diagrams in one-loop MHV 
amplitude. Of course, the one-loop MHV amplitude is known analytically, but this result 
will be extremely useful in analyzing the diagrams contributing to the two-loop MHV 
amplitude. 



A. 1.3 Differential equations for the one-loop MHV amplitude 

So far, we analyzed the one-loop amplitude giving the equations for all diagrams contribut- 
ing to it. However, there exists an equation for the integrand of the whole amplitude. In 
fact, for the n-point amplitude there exist n 2 differential equations of the form 

Qij^MHv = i,j = l,...n (A. 165) 

where 

n 

Qn = 0ik0k i - ( n - 2 ) 0i i ■ (A.166) 
fc=i 

Note that the structure of operators is quite similar to that which we found for the pentagons 
at the end of the last subsection. 



A. 2 Diagrams in the two-loop amplitude 

We start with the form of two-loop MHV amplitude given in [T3] and reviewed in section 
[2J It contains three topologies: double-boxes, penta-boxes and double-pentagons. 

We will analyze each topology separately. Naively, they look like one-loop diagrams 
(two-mass hard integrals and pentagons) glued together, and we will see that in some sense 
they really are. 

We want to find a second-order differential equation for the integrands of these diagrams, 
keeping in mind our motivation that acting with the correct second-order differential op- 
erator on a two-loop diagram can decrease the loop level by one). Therefore, it would be 
very natural to act just on a one- loop sub-diagram, which is exactly what we are going to 
do. We will see that there is close link with the operators found in the last subsection. 

A. 2.1 Double-boxes 

The integrand for the general double-box appearing in the amplitude 

ji _ (1234)(z-l^+H+2)(23^+l) 

double-box (AB12)(AB23)(ABM)(ABCD)(CDi-li)(CDii+l)(CDi+li+2)' K ' ' 



can be rewritten in a very suggestive form 

(1234)(23 i i+l) ii+1 i+2)(AB 



(AB12) (AB23) (AB3A) (AB (ABCD) (CD i-1 i) (CD i i+l) (CD i+l i+2) 



. (A.168) 



We see that the double box can be written in a factorised form as a product of two-mass 
hard diagrams (note that they are properly normalised), where the right CD sub-diagram 
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contains twistors A,B as external, but the left AB sub-diagram depends only on external 
twistors. If5<i<n— 1 (both internal corners are massive), then there is no overlap 
between AB and CD sub-diagrams and the AB part is annihilated by the same operators 
as the normal 2-mass-hard integrand, 

Ol224,hbox = O^ihbox = Oi23lhbox = O '433lhbox = (A. 169) 

where here by the normalisation iV we mean the normalisation of left sub-diagram N = 
(1234) (23 i z+1) . If i — 4, 5, n — 1, n just half of the operators survive (we still have one 
massive corner). 

For n = 4, 5, 6 we have diagrams that do not contain any massive corners. However, 
the diagrams for n = 5 and for n = 6,i = 5 satisfy differential equations and here we use 
the complicated forms that we found in the end of the discussion of the two-mass hard 
boxes. Thanks to restrictions we made on the first-order operators for the construction of 
the second-order operator, these operators act just on one half of the double-box integrand 
even if the indices overlap (which exactly happens for i = 5 or i = n — 1). The only 
constraint is that at least one external line (which can be massless and just on one side 
of the diagram) must be attached at the point were the two one-loop diagrams are glued 
together. Then the equations are exactly the same as we found for the two-mass hard case, 
since we act on just half of the diagram. 



A. 2. 2 Pentagon-boxes 

As in the previous case, we can write the integrand for the pentagon-box appearing in the 
two-loop MHV amplitude, 

jpenta-box = (AB(123) H (i-li i +1) ) (2 i J J +1 ) (j-l j j +1 j +2) 

n;i ' j (ABl2)(AB23)(ABi-li)(ABii+l)(ABCD)(CDj-lj)(CDjj+l)(CDj+lj+2) ,[ ' ' 

in a factorised form, 

(45(123) n (i-lii+l))(2ij j+1) (j-ljj+lj+2)(ABjj+l) 



(ABU) (AB23) {AB i-1 i) {AB i i +1) {AB j j +1) (ABCD) (CD j-l j) (CD j j +l)(CD j +1 j +2) 

(AB(123)n(i-lii+l))(CD2j) (2 z j j +1) (j-l j j +1 j +2) 

(ABl2)(AB23)(ABi-li)(ABii+l)(ABCD) ' (CD2j)(CDj-lj)(CDjj+l)(CDj+lj+2) ' 



(A.171) 



We see that there are two choices, either the pentagon or the two-mass-hard is the "proper" 
one- loop sub diagram. 

If we use the first choice, then the the two-loop diagram is annihilated by all operators 
that annihilate its pentagon one-loop sub-diagram, and the same is true for the second 
choice, where we use the differential operators for the two-mass-hard diagram. If the indices 
do not overlap, we can just act with the operators that we have for pentagons, resp. two- 
mass hard and we annihilate the integrand. In the special case of j = i + 1 or j = n — 1 
we would repeat the same discussion as in the case of double-boxes. Here we use the more 
complicated forms that we derived in the end of the two-mass hard and pentagon sections 
that can be used even if the indices overlap. Therefore, all pentagon-box integrands satisfy 
second-order differential equations. 
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A. 2. 3 Double-pentagons 

Finally, we can write the integrand for the double-pentagons, 

(AB(i23)n(£-iee+i))(CD(i-iu+i)nti-ijj+i))(2£ij) 

(ABl2)(AB23)(AB£-l£)(AB££+l)(ABCD)(CDi-li)(CDii+l)(CDj-lj){CDjj+l)' [ ' ' 

in a factorised form, 

(AB(123) n (£-1 ££+!)) (2 £i j) (CD(i-lii+l) n jj+l)){ABij) 

(AB12) (AB23) (AB £-1 £) (AB £ £+l){AB i j) (ABCD) (CDi-li)(CDii+l)(CD j-1 j)(CD j j +1) ' 

We again see two correctly normalised pentagons, and one of them is a proper pentagon on 
which we can apply differential operators. The discussion would be identical to the previous 
cases, one just uses the operators found in the pentagon section which are guaranteed to 
annihilate also the double-pentagon integrand. 

B Five-point pentagon and penta-box integrals 

As an explicit example of the differential equations for regulated integrals of section [3J, we 
give the result for the five-point pentagon and penta-box integrals at one and two loops. 
Let us start with the one-loop pentagon integral. Without loss of generality, we choose it 
in the orientation 

t [ d*Z AB (2345) (4512) (AB13) 

5;5 ' 1,3 J m 2 (AB51)(AB12)(AB23)(AB34)(AB45) ' 

and we recall the reader that because of infrared divergences we also include +m 2 terms 
in the definition of the propagators, see section [2j Introducing Feynman parameters, we 
obtain 

F^ = ~24A4, f°° +OK). (B.174) 

The small m 2 expansion of this integral can be conveniently obtained using Mellin-Barnes 
methods. The answer depends on three variables, yi = xl 5 /xl 3 ,y 2 = x\^jx\^y 3 = x1 4 /m 2 , 
and can be written as 

FlS,s = -\ log 2 (2/iy 2 2/ 3 ) - 2Li 2 (l - yi) - 2Li 2 (l - y 2 ) + ^ + 0(m 2 ) . (B.175) 
2 o 

We can also compute the corresponding two-loop penta-box integral. In the nomenclature 
of section [21 it is 

F penta-box = f d'Z AB d'Z CD (m 2 ) ~ 2 N(CD 13) 

5;5,i,3 J (AB34){AB45){AB5l){ABCD){CD5l){CD12){CD23){CD34:) ' 1 ' ' 
where N = (3451) (2345) (4512). It evaluates to 

^Kr b ° X = Kvmys) ~ 2H (y iym ) {f( yi ) + f(y 2 )) + \g{ Vl ) + g(y 2 )) + 0(m 2 ) ,(B.177) 
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where 

h(x) = -^ i + 2C,H {x)- l -'K 2 H 0fl {x)-H Qfiflfi {x), (B.178) 

g(x) = ~^H 2 (x) + 2H (x)H 2fi (x)-AH 3fi (x)-2H 2A0 (x)-AH 2Afi (x), (B.179) 
f{x) = #o,i,o(z). (B.180) 
The integrals satisfy 

y^y^/^l = - 1 + . (B.i8i) 

and 

l/i^l/2^^5t3" bOX = *SSs + 0(m 2 ) . (B.182) 
One can check that equations ( 13.971) and ( 13.99j) are indeed satisfied. 
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